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Abstract. The internal structure of a slowly rotating, charged black hole that is 
undergoing mass inflation at its inner horizon is derived. The equations governing 
the angular behavior decouple from the radial behavior, so all conclusions 
regarding inflation in a spherical charged black hole carry through unchanged for 
a slowly-rotating black hole. Quantities inflate only in the radial direction, not in 
the angular direction. Exact self-similar solutions are obtained. For sufficiently 
small accretion rates, the instantaneous angular motion of the accretion flow has 
negligible effect on the angular spacetime structure of the black hole, even if the 
instantaneous angular momentum of the accretion flow is large and arbitrarily 
oriented. 
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1. Introduction 

The internal structure of a realistic (astronomical) rotating black hole has remained 
an unsolved problem for decades (e.g. §20.5 of jl], §5.3 of [2]). The structure of a 
rotating black hole in the process of collapsing is undoubtedly thoroughly complicated 
[3]. However, thanks to the no-hair theorem [H [3 [6], it is expected that a newly 
collapsed rotating black hole will quickly relax to the Kerr geometry. Aside from the 
initial collapse event, and rare occasions of high accretion such as a black hole merger, 
a real astronomical black hole spends the great majority of its life accurately described 
by the Kerr geometry down to near its inner horizon. 

The reason that the Kerr geometry fails inside a black hole is that at (or rather 
just above) the inner horizon, the Kerr geometry is subject to the relativistic counter- 
streaming instability discovered by [7j, and called by them mass inflation. The 
inflationary instability is the nonlinear realization of the infinite blueshift at the inner 
horizon first pointed out by [8]. Mass inflation inside spherical charged black holes 
has been investigated extensively (see e.g. the lists of references in ^ US]), but there 
have been few studies of inflation inside rotating black holes [lOl [Til HH US] (see also 
the review [14]). 

The present paper addresses, and solves, the problem of the internal structure of 
a slowly rotating, charged black hole that is undergoing inflation at its inner horizon. 
Inflation requires the simultaneous presence of ingoing and outgoing streams near the 
inner horizon; but even the tiniest incident streams suffice to induce inflation. The 
paper follows earlier papers in this series [15[ I16[ 117] in assuming that ingoing and 
outgoing streams are being generated continuously and steadily by accretion. This 
assumption differs from most of the literature on inflation, which has focused on the 
situation where a black hole collapses and then remains isolated [18j . in which case 
the outgoing stream is dominated by a decaying Price tail [121 HH] of gravitational 
radiation. Real astronomical black holes are however never isolated, and it seems 
likely that any decaying Price tail produced by the initial collapse event will soon be 
overwhelmed by streams generated by accretion. Matter that falls from outside the 
horizon is necessarily initially ingoing, but it may become outgoing inside the horizon 
if it has enough charge (if the black hole is charged) or enough angular momentum (if 
the black hole is rotating). 

The approach followed in the present paper is to introduce a small rotational 
perturbation to a spherically symmetric black hole that is undergoing inflation. It is 
necessary to assume that the black hole is charged so that the unperturbed geometry 
has an inner horizon where inflation is taking place. In this sense the solution is not 
realistic, since a real black hole probably has tiny charge. However, one may hope 
that the solution for the internal structure of a slowly rotating, charged black hole 
may shed light on the real problem of a fully rotating, uncharged black hole. 

As shown by TfJ, the radial length scale over which the interior mass and other 
gauge-invariant quantities (such as the Weyl curvature, and the center-of-mass density 
and pressure) exponentiate during inflation is proportional to the accretion rate of 
the incident ingoing and outgoing streams. Thus, counter-intuitively, the smaller 
the accretion rate, the more rapidly inflation exponentiates as a function of radius. 
In the limit of a tiny accretion rate, the structure resembles a step function at 
the inner horizon, above which the Kerr-Newman geometry applies to an excellent 
approximation, and at which the interior mass and other quantities inflate to absurdly 
huge values while the radius scarcely changes. 
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Astronomically realistic black holes have small accretion rates over the great 
majority of their lifetimes, the accretion timescale being much longer than the black 
hole crossing time. For small accretion rates it might seem justifiable to assume 
the stationary approximation d/dt = (which [551 [23j terms the homogeneous 
approximation, since the Killing vector associated with time translation symmetry 
is spacelike inside the horizon). It is true that the spacetime is almost stationary, 
temporal partial derivatives being much smaller than radial partial derivatives, d/dt <^ 
d/dr, but it is strictly stationary only in the unrealistic case of symmetrically equal 
ingoing and outgoing streams [17j . The reason that stationarity breaks down is that 
quantities such as the interior mass M inflate so rapidly that dM/dt is never small in 
the inflationary zone even in the limit of infinitesimal accretion rates. 

Although the slowly-rotating black hole spacetimes considered in this paper are 
not time-translation invariant, the equations do admit solutions with conformal time- 
translation symmetry, or self-similarity. Self-similar solutions are considered herein. 



2. Slowly rotating equations 



2.1. Slowly rotating perturbation 

This paper follows [HI [TH [17] in adopting the following line-element for a general 
spherically symmetric background spacetime, in polar coordinates = {t, r, 9, 0}: 

ds2 = -^^ + -^ fdr-A-V+r2(d02+sin20d<^2) 
\ a J 

Through the identity 

ds^ = g,,, Ax^'dx" = we^^^e", Ax'^Ax'^ , (2) 

the line-element ([IJ encodes not only a metric g^^, but a complete inverse vierbein 
e™^, corresponding vierbein em'', and orthonormal tetrad {7t, 7^., 7^, 7^}, satisfying 
7,„ • 7„ = rjj-an with ?7„i„ the Minkowski metric. The coefficients in the line- 
element ([!]) constitute the components of a covariant tetrad-frame 4- vector, the radial 
4-gradient 

/3m = dynT , (3) 

where dm = 7m • d = e„i^d/dx^^ denotes the directed derivative along the 7^ 
tetrad axis (not to be confused with the coordinate derivative d/dx^). Physically, 
the directed derivatives dt and dr are the proper time and radial derivatives measured 
by a person at rest in the tetrad frame. The scalar product of j3m with itself defines 
the interior, or Misner-Sharp [24j . mass M, a gauge-invariant scalar 

^-l^-/3m/3™=/3,2-/?,2 . (4) 

Details of the vierbein, tetrad-frame connections, and tetrad-frame Riemann, Einstein, 
and Weyl tensors that follow from the unperturbed line-element ((T|) are summarized 
in Appendix A of [17 . 

Perturbations to the spherically symmetric spacetime can be characterized by the 
dimensionless covariant tetrad-frame vierbein perturbation ^pmn defined by 

e™^ = (C - ^m")en'' , (5) 
where the overscript " signifies an unperturbed quantity. Because the perturbation 
^Pmn is of linear order, it is legitimate to raise and lower its indices with the 



The interior structure of slowly rotating black holes 



4 



unperturbed tetrad-frame metric (the Minkowski metric), and to transform between 
coordinate and tetrad frames with the unperturbed vierbein e„i^ and its inverse e™^. 
Thus the vierbein perturbation ipmn can be regarded as a tensor field defined on the 
unperturbed sphericaUy symmetric background. 

Spherically symmetric spacetimes are invariant under the 3-dimensional group 
0(3) of orthogonal transformations. It is natural to expand perturbations in 
eigenmodes of that group. A perturbation eigenmode is characterized not only by its 
harmonic numbers Im, but also by its spin s, which expresses how the perturbation 
transforms under rotations about the radial direction. Perturbations ipmn admit scalar 
(spin-0), vector (spin-±l), and tensor (spin-±2) modes. The spin of a perturbation 
ipmn is most easily identified not in an orthonormal tetrad, as here, but rather in a 
spinor tetrad, where the angular tetrad axes 'jg and 7^ are replaced by spinor axes 

7+ = -^(79 + 17^) , 7- = -^(70 -170) • (6) 

The general rule is that the spin of a tensor equals the sum of the +'s and — 's of its 
covariant indices in the spinor tetrad. 

A slowly rotating spacetime is obtained from a spherically symmetric spacetime 
by a perturbation that is (a) cylindrically symmetric, (b) axial, meaning that the 
perturbation changes sign under a flip 7^ — > — 7^ of the azimuthal tetrad axis, and 
(c) dipole, the lowest order non-spherical harmonic. The slowly rotating conditions 
imply that the vierbein perturbation ipmn is a spin-1 dipole harmonic, 

(Pmn OC iFio(0, 0) = . (7) 

VOTT 



Of the 16 components of the 4x4 tensor of perturbations ipmn, there are 4 that satisfy 
the conditions of being spin-1 and axial, namely 

Vt0 I 'Pr4> , Vcfit , fcfir ■ (8) 

The spin-1 axial character of these perturbations is manifest in a spinor tetrad, where 
the possible spin-1 perturbations are the 4 complex quantities (^t+, fr+, f+t-, and f+r, 
and the axial condition implies that each of these 4 perturbations is purely imaginary. 

Of the 4 possible degrees of freedom represented by the perturbations (O, only one 
degree of freedom is coordinate and tetrad gauge- invariant, and therefore represents 
a real physical perturbation. The remaining 3 degrees of freedom represent gauge 
freedoms, comprising an infinitesimal coordinate transformation of the azimuthal 
coordinate (f>, and infinitesimal tetrad (Lorentz) transformations in the t-(j) and r-(j) 
planes. 

Under an infinitesimal coordinate gauge transformation of the azimuthal 
coordinate (p, 

cl)^(f> + e, (9) 
the vierbein perturbations transform as 

/ de de\ 

ft4> ^ <ft,p-r sm 9 la— + Ptg;^] , (10a) 
de 

iprcji ^ frcji - r sin 9 Pr-g;^ , (106) 
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with if^t and ipif,r being coordinate gauge- invariant. If the spacetime were either time 
translation-invariant or conformally time translation-invariant, so that de/dt = 0, 
then the combination 



(11) 



would be coordinate gauge-invariant. In general the spacetime is not time translation- 
invariant, conformal or otherwise, but as commented in the Introduction it is 
legitimate to think of the spacetime as being almost stationary, and therefore the 
combination (jlip is almost coordinate gauge-invariant. 

The antisymmetric combinations Lpttp — ^ci>t and (pr^ — ip^r are not tetrad gauge- 
invariant, the former transforming under an infinitesimal Lorentz boost in the t-4> 
plane, the latter transforming under an infinitesimal spatial rotation in the r-cj) plane. 
The symmetric combinations ipt^, + (pct,t and iprcj, + ^4,7- are tetrad gauge- invariant. If 
the spacetime were time translation-invariant or conformally time translation-invariant 
then the combination 

Prift^, + f^t) - fiti'Pr,!} + 'P4>r) (12) 

would be gauge-invariant with respect to both infinitesimal coordinate and tetrad 
gauge transformations. Again, in general the spacetime is not time translation- 
invariant, but it is legitimate to think of the combination ([T^ as being almost 
coordinate and tetrad gauge-invariant. In terms of the almost-gauge-invariant 
perturbation oj defined below, the combination (|12p equals — ra/3rW sin^^. In the 
general case, there is no linear combination of vierbein perturbations ipmn that is 
fully gauge- invariant, but the axial spin-0 component of the Weyl tensor provides 
a differential combination W of vierbein perturbations, equation ([22|l . that is fully 
gauge- invariant . 

In view of the behavior under gauge transformations just discussed, it is 
convenient to introduce 4 quantities to, ip, C*j and C representing the different types 
of perturbation as follows: 

Lo almost gauge-invariant , 
ip coordinate transformation of , 
Lorentz boost in the t-(j) plane , 
spatial rotation in the r-(f> plane , 
in terms of which the vierbein perturbation tpmn is 

/ -raw-/3tV + C' \ 



(13a) 
(136) 
(13c) 
(13rf) 



sin 9 











(14) 



The perturbations are tetrad-frame components of the Killing vector associated 
with cylindrical symmetry, equation (jlSp . The resulting perturbed vierbein is 



a 




f3t 




y -C*Q!sin6' (3^ 




1 

rsin^ 



(15) 
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and the perturbed inverse vierbein is 



1 

a 









rC* sin^ 


1 

Wr 






r 


/■r ■ 2 

rQ sm 



rfl sin 9 


-5' sin 61 





rsinO 



(16) 



where 



f2 = w - 



raPr 



* = ■(/) + 



(17a) 
(176) 



The final column of the inverse vierbein ([TC]) constitutes the spacelike 
contravariant tetrad- frame Killing vector rsin^^C™ associated with cylindrical 
symmetry 



^^e^.d^=rsm^erd„ 



(18) 



The second column of the vierbein constitutes the covariant tetrad-frame 
4- vector (3rm the radial 4- gradient defined by equation ([3]). The azimuthal component 
Pci> of the 4-vector follows from dr/d(j> = 0, or equivalently C™/^m = 0, which implies 

f3^ ^~{Cf3t + Cf3r) sine . (19) 

The inverse vierbein (|16p is equivalent to the line-element 

2 / , , n \ 2 



ds^ 



di 



q2 I „;„2 



+ C* sin' 



+ 



dlnr 



-I- C'^ sm ^ 

Pr ra 



de^ + sin^e (d<?!) - * dlnr - r2 dt)' 



(20) 



2.2. Einstein and Weyl tensors 



From the line-element ([201) can be derived the tetrad- frame connections, and tetrad- 
frame Riemann, Ricci, Einstein, and Weyl tensors in the usual way. In the present 
case of a spherical spacetime subjected to an axial vector perturbation, all quantities 
can be classified as either scalar (spin-0) or vector (spin-1). 

The only spin-0 component of the Riemann tensor that differs from its value in 
the unperturbed spherical background is the axial part of the spin-0 component of 
the Weyl tensor C'kimn ■ The component is proportional to the spin-0 dipole harmonic 
0^10 = cos 6'/-\/47r, and is conveniently written 



Wcos9_ 1 

5 = ^-Z^tred) — ^^tBr<h ~ ^tdirB 

2 



(21) 



where W{t,r) is a function only of t and r, not of angular variables. In terms of the 
vierbein perturbations, W is explicitly 

2 Vainr ~ 'dt 

The axial spin-0 component of the Weyl tensor, equation (I^T]) . hence W, is fully gauge- 
invariant with respect to all relevant gauge transformations, which is to say not only 



W 



(22) 
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infinitesimal coordinate and tetrad transformations, but also arbitrary radial Lorentz 
boosts. 

There are 4 non- vanishing axial spin-1 components of the tetrad- frame Riemann 
tensor Rkimn, yielding the following 4 non- vanishing spin-1 components of the tetrad- 
frame Einstein tensor Gmn and Weyl tensor Ckimn- 

Gt0 = — {Rtee<j> + Rtrr4>) , (23a) 

Gr4, = — {Rr004> + Rtrt<t>) i (236) 

Ctee4> = —Ctrr<i> = ^ {Rtee<t> — Rtrr<j>) , (23c) 



Cr 



Gtrt<t> — o {RrOe4> — Rtrt4 



(23 d) 



Being tetrad-frame quantities, all these components are automatically coordinate 
gauge-invariant. However, none of them arc tetrad gauge- invariant: all spin- 
1 components of the Riemann tensor transform under infinitesimal spin-1 tetrad 
(Lorentz) transformations. 

The non- vanishing spin-1 components of the tetrad-frame Einstein tensor Gmn 

are 



Gt4, = sin 6 
Grd, = sin t 



^drir^W) - C{Gu + G^^) - CGt 



\dt{r''W)-C{Grr- 



- CGt: 



(24a) 
(246) 



Note that the derivatives dm m equations (24) are directed (tetrad-frame) derivatives, 
not partial (coordinate-frame) derivatives. 

The non- vanishing spin-1 components of the tetrad- frame Weyl tensor Gkimn are 



Gtee4, = sin 6 
Cr666 = sin ( 



drW - 3C*C 
dtW + 3CC 



(25 a) 



(256) 



where C is the Weyl scalar, the polar part of the spin-0 component of the Weyl tensor 

-Gt<j,t<p 



C: 



trtr 



teie 



-CrOrB — Cfihrd 



g (Gtt — Grr + 



M 



(26) 



The Weyl scalar C is invariant with respect to arbitrary radial Lorentz boosts. 



2.3. Electromagnetic field 

The black hole is necessarily charged, in order that the unperturbed black hole have 
an inner horizon where inflation can take place. 

The unperturbed spherical background admits just one component of the 
electromagnetic field Fmn, the polar spin-0 component, the radial electric field 
Er = Ftr. This is a monopole field, and it can be written 



Er = 



(27) 



which defines the charge Q{t, r) interior to radius r. 
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The rotational perturbation introduces scalar (spin-0) and vector (spin-1) axial 
dipole components to the electromagnetic field. The axial spin-0 part is the radial 
magnetic field B^. = F^g . The component is proportional to the spin-0 dipole harmonic 
oYio = cos9/\/iTT, and can be written 

B cos 6* 

Br ^ , (28) 

which defines Bit, r). Both Q and B, equations ([27|) and ([28]) . are fully gauge-invariant 
with respect to infinitesimal coordinate and tetrad gauge transformations, and to 
arbitrary radial Lorentz boosts. 

The non- vanishing axial spin-1 components of the electromagnetic field are 
Etj, = Ft^ and Bg = Fr^. The source-free Maxwell's equations £^^"^^DiFmn = 
imply that the axial spin-1 components of the electromagnetic field must be related 
to Q and B by 

2r 



E^^-s\ne[—dtB + C^^], (29a) 



B,^_sin(?(la,.i?-C*§) . (29&) 



The source Maxwell's equations D^'^Fmn = 47rj„ relating the electromagnetic field 
Fmn to the electric current jVi are 

d,Q = ^inr^jt , (30a) 

dtQ = -Wj, , (30&) 

i (^-D'dt - D^-dr + ^) S + = 47rrC"j™ , (30c) 

where Dm signifies the tetrad-frame covariant derivative. 

The electromagnetic energy-momentum tensor is given by 47rT^„j = FmkFn^ — 
jrjmnFkiF^^ . In terms of the electric and magnetic fields (|27p and (29), the axial 
spin-1 components of the electromagnetic energy- momentum tensor are 

AnT^^ - -ErBe , 4^r;^ = -E^E^ . (31) 

Explicitly, from equations (29) inserted into equations ((3T|) . 

87rT/^ = sin0(^^a,B + C*^^ , (32a) 
SttT,^. = sin(?f%B-C''^) . (32&) 



The covariant divergence of the electromagnetic energy-momentum is Z),„TJ"" = 
-F"™jm, and the component of this along the azimuthal Killing vector is 

CuDmTr = ^rdmB . (33) 
Ir 

2.4- Perfect fluid streams 

The inflation instability requires the simultaneous presence of ingoing and outgoing 
streams near the inner horizon [7] . I follow |16lll7j in considering two separate streams, 
charged "baryons" with a relativistic equation of state (a value slightly less than 1/3 



The interior structure of slowly rotating black holes 



9 



allows for the expected slow increase with temperature of the effective number of 
relativistic species), 

w = Pb/pb = 0.32 (baryons) , (34) 

and neutral "dark matter" with a pressureless equation of state, 

Pdl Pd — (dark matter) . (35) 

Repelled by the charge of the black hole generated self-consistently by their accretion, 
the baryons become outgoing inside the horizon. The dark matter falls freely into the 
black hole, and remains ingoing. Inflation is produced by relativistic counter-streaming 
between the outgoing baryons and the ingoing dark matter. 

Each of the two streams constitutes a perfect fluid. The equations in this 
subsection are for a general perfect fluid. 

The tetrad-frame energy-momentum T™" of a perfect fluid of proper density and 
pressure p and p moving at tetrad-frame 4-velocity u™ is 

T"'" ^{p + p) u"'u'' + p r/™" . (36) 

Since the spacetime is cylindrically symmetric, any geodesic is characterized by 
a conserved azimuthal angular momentum, a gauge-invariant scalar. The azimuthal 
angular momentum (per unit mass) of a stream is given by the covariant azimuthal 
component of the coordinate-frame 4-velocity (the coordinate-frame 4-velocity 
t!^ is written with a Greek upsilon to avoid any possible confusion with the tetrad- 
frame 4-velocity u"\ written with a Latin u): 

= e"^ nmn = r sin^ Cnu"' ■ (37) 

It is convenient to rewrite the azimuthal angular momentum of a stream in terms 
of an angular momentum parameter a defined by 

U0 = asin^^ . (38) 

If the stream is slowly rotating, and if it has no 9 motion (it^ = 0), so that its rest 
frame belongs to the system of slowly-rotating tetrads considered in then a = r^t 
in the tetrad rest-frame of the stream. However, equations ([57)) and (|38p remain valid 
in general, regardless of the angular motion of the stream. Equations (|37p and ([55)1 
imply that the azimuthal component u*^ of the tetrad-frame 4-velocity of a stream is 
related to the time and radial components m* and u^, and to its angular momentum 
parameter a, by 

u'* = sin6l(--Ctu*-Crw'') . (39) 

With u'^ related to u* and by equation ([5^ . it follows from equation ([55)1 that 
the axial spin-1 components of the tetrad- frame energy-momentum tensor of a perfect 
fluid are 

a{p + p)u^ 



rpt4> _ ^ 



smf 



rpr4> 



r 



(40 a) 



(406) 



r 

The component of the covariant divergence of the perfect-fluid energy-momentum 
along the azimuthal Killing vector C" is 

C„^™r™" = 1 { (a* -f hr) [r^a{p + py] + Or + ht) [r^a{p + p^] } , (41) 
where km = ^trm are tetrad- frame connection coefficients, equation (64) of |17j . 
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2.5. Einstein equations 
The Einstein equations are 

Gmn = SnTmn ■ (42) 

With the Einstein tensor from equation (24) and the electromagnetic and perfect fluid 
energy-momentum tensors from equations (|3ip and (40), the axial spin-1 Einstein 
equations reduce to 

drir^W) ^ QdrB - Snr^ ^ a{p + p)u^ , (43a) 

streams 

dt{r^W) ^ QdtB + Sirr^ J2 a{p + p)u'' , (43&) 

streams 

where the sum is over the accreting baryonic and dark-matter streams. These 
equations (43) are gauge-invariant, the parts depending on the gauge perturbations 
and in equations (24), ([?!]) . and (40) having canceled between the left and right 
hand sides. 

If the underlying sources of energy-momentum T,„„ are arranged to satisfy 
conservation of energy-momentum, as they should, then one of the two Einstein 
equations (43) can be treated as redundant, since it serves simply to enforce overall 
covariant energy-momentum conservation in the azimuthal direction 

rO^Gran - . (44) 

That there is effectively only one independent axial spin-1 Einstein equation is 
consistent with the fact that there is only one physical degree of freedom in the axial 
spin-1 vierbein perturbations. 

2.6. Energy exchange 

To complete the equations, it is necessary to specify the exchange of energy between 
each pair of species. 

I assume for simplicity that the dark matter and the baryons stream freely through 
each other without interacting. 

For the interaction between charged baryons and the electromagnetic field, I 
assume Ohm's law in the rest frame of the baryons, 

jm = crEm (m = r, 0) , (45) 

where a is an electrical conductivity. The azimuthal current r^^jVri that appears on 
the right hand side of the Maxwell equation (|30cp is then, in the rest frame of the 
baryons, 

rrjm = -[ajt + ^dtB^ . (46) 
3. Slowly-rotating Kerr-Newman 

As long as the accretion rate is small, as is true for astronomically realistic black 
holes, the spacetime down to near the inner horizon is well-approximated by the 
slowly-rotating Kerr-Newman [551 US] geometry. The Kerr-Newman geometry is 
characterized by the mass M,, charge Q,, and angular momentum parameter a, of 
the black hole, all constants. 
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The monopole part of the slowly-rotating Kerr-Newman solution coincides with 
Reissner-Nordstrom. The interior mass M(r) and charge Q{r) are 

M = M,-^ , Q = Q, . (47) 
2r 

If the time coordinate t is chosen to reflect the stationarity of the spacetime, d/dt — 0, 
and scaled globally to synchronize with proper time at rest at infinity, then the vierbein 
coefficient a is related to the vierbein coefficients f3m by 

The radial 4-gradient f3m, equation ([3]), is a tetrad-frame 4- vector, and different choices 
of f3t or Pr correspond to different gauge choices of radial Lorentz boost of the tetrad 
frame. For example, Pt — recovers the standard diagonal Reissner-Nordstrom line- 
element, while /?r = 1 yields the GuUstrand-Painleve [27] line-element. 

The dipole part of the slowly-rotating Kerr-Newman geometry is specified by 
the gauge-invariant vierbein perturbation Lu{r), and by the gauge- invariant quantities 
W{r) and B{r) describing the axial spin-0 Weyl and electromagnetic scalars, 
equations and ((28)) : 

oj = — ^ , (49a) 

W^=^(-3Af. + ^) , (49&) 

B = ^ . (49c) 
r 

The Kerr-Newman geometry is most commonly expressed in the Boyer-Linquist |28j 
tetrad, where the spin-1 components of the Einstein, Weyl, and electromagnetic 
tensors all vanish identically. The Boyer-Linquist tetrad is attained with the gauge 
choices 

rC = ~a,Pr , rC = a,Pt ■ (50) 
The coordinate gauge perturbation i/;, or equivalently ^, equation (jl76|) . plays no role, 



but it simplifies the line-element ((20)) to set 

* = . (51) 



4. Self-similar solutions 



The slowly- rotating black hole equations admit self-similar solutions. Though idealized 
as a model of reality, the similarity solutions are exact solutions that provide insight 
into the general case. The similarity solutions in this section are extensions of those 
discussed by [HldSKn]. 

As has been seen in ^ the equations governing the slowly-rotating perturbation 
are completely decoupled from the spin-0 monopole equations governing the spherically 
symmetric background. It follows that the behavior of inflation in slowly-rotating 
black holes is identical to that in spherical black holes. This conclusion is true 
subject to the consistency condition that the slowly-rotating perturbation remains 
small through inflation. As seen in i j4.4l below, the slowly-rotating perturbation 
remains small, so the solutions are indeed self-consistent. 
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The boundary conditions of the sphericaUy symmetric self-similar solutions are set 
at a sonic point outside the outer horizon, where the infalling baryonic fluid accelerates 
from subsonic to supersonic. The mass M, and charge Q, of the black hole can be 
defined to be the mass and charge that a distant observer would see if there were no 
matter or charge outside the sonic point: 

n2 

Af, = M + , Q, =Q at the sonic point . (52) 
2r 

The extra mass /2r added to the interior mass M is the mass-energy in the electric 
field outside the sonic point, given no charge outside the sonic point. 

For the spherically symmetric solutions, there are three dimensionless boundary 
conditions at the sonic point: (1) the mass accretion rate M,; (2) the charge-to-mass 
ratio Qt/Mt] (3) the dark-matter-to-baryon density ratio pd/pb- The solutions scale 
with the mass M, of the black hole. 



Similarity assumption 

Self-similarity is the assumption that the system possesses conformal time translation 
invariance, that is, there exists a conformal time t such that the system at any one 
conformal time is a scaled copy of the system at any other conformal time. 

The line-element (|20l) can be cast in a form that has explicit conformal symmetry 
by transforming the time and radial coordinates t and r to conformal time and radial 
coordinates t and defined by 

t = T, lnr==_R + T. (53) 

The fact that it is possible to choose the time coordinate to be conformal time, t — t, 
reflects the gauge freedom in the choice of time coordinate. In terms of the conformal 
coordinates t and R, the line-element pO|) is 

2 



ds^=r^- {edT + Csm^0dcf>Y +(^+CdT + CsinH 



' dB. 

Jr 

+ de^ + [sm 9 {d<t) ~ ^ dR) dr]^^ . (54) 

The circumferential radius r (squared) appears as an overall conformal factor. 
Conformal time translation symmetry requires that all conformal inverse-vierbein 
coefhcients, the coefficients of differentials inside the braces in equation ([M]) . are 
functions only of conformal radius fl, not of conformal time t. 

The in the line-element ([Si]) are the tetrad-frame components of the Killing 
vector associated with conformal time translation invariance, 

|- = e", dm = rrdm . (55) 

In terms of the vierbein coefhcients of the line-element (|20|) . the conformal Killing 
vector ^"^ is 

r^{;^,^(-^ + l),0,-(^^ + *)sm.}. (56) 

The scalar product of the conformal Killing vector with itself defines the horizon 
function A, a coordinate and tetrad gauge-invariant scalar, 

A ^ = e - • (57) 
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The horizon function is positive outside the horizon, zero at the horizon, and negative 
inside the horizon, reflecting the fact that the Kilhng vector ^™ is respectively timeUke, 
null, and spacelike outside, at, and inside the horizon. The scalar product of the 
conformal Killing vector ^™ with the radial 4-gradient f3rn follows from dlnr/dr — 1 , 
which implies 

rP,n = 1 . (58) 

In self-similar solutions, all quantities are proportional to some power of 
circumferential radius r, and that power can be determined by dimensional analysis. 
Dimensional analysis shows that the conformal coordinates {t, _r, 9, cj)}, the coordinate 
metric (7^,^, and the vierbein perturbation (fmn £^re all dimensionless. In particular, 

(3„i , S,"^ , uj , ip , are dimensionless . (59) 

The vierbein e„i'^, inverse vierbein e™^, tetrad- frame connections Tkmm and tetrad- 
frame Riemann tensor Rkimn, energy- momenta Tmn, and electromagnetic field Fm„, 
scale as 

OC r , 6 OC r , ^kmn CX] r , 
Rklmn OC r"^ , Tmn OC r~^ , Fmn OC r~^ . (60) 

Various specific quantities that appear in HA.'I] scale as 

W (X r'^ , Q OC r , B oc r , a oc r . (61) 



4-2. Similarity equations 

Four equations ([M)) - (p7|) govern the angular structure of the slowly-rotating black 
hole, one for each of the baryons, dark matter, gravity, and electromagnetic field. 
For simplicity, the equations in this subsection are for zero conductivity, cr = 0, as 
assumed in the model of ii4.4l Equations for finite conductivity are given for reference 
in Appendix I Appendix A[ 

I follow [15l [16] in choosing the integration variable x to be the dimensionless 
baryonic time defined by da; = r^^dr, where r is the proper time along the worldline 
of the baryons. Equivalcntly 

^ = rdt (62) 

ax 

where dt is the directed time derivative in the rest frame of the baryons. The 
dimensionless baryonic time x increases monotonically, since proper time does. The 
dark matter falls into the black along a different trajectory than the baryons, and 
it is therefore necessary to distinguish the dark matter radius and dark matter 
tetrad frame from the baryonic radius and tetrad frame, as detailed in '16]. I treat 
the baryonic frame as the default frame, and for brevity drop the baryonic subscript 
b from r and ^™ evaluated in the baryon frame. 

In self-similar solutions, the volume of a Lagrangian volume element is 
proportional to ^Trr^i^''. Angular momenta L of Lagrangian volume elements can 
be defined to be L = ^Trr'^£^^a{p+p). Explicitly, for each of the baryon (subscripted b) 
and dark matter (subscripted d) streams, whose equations of state are given by 
equations and 

Lb = ^nr^Co-bi^ + w)pb , (63a) 
Ld = |7rr^^S«dPd . (636) 
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Covariant conservation Cn-Dm(T"™ + T™") = of electromagnetic plus baryon energy 
along the azimuthal Killing vector C™, equations ([55]) and (|^T|l . yields the following 
equation of conservation of angular momentum of baryons in their rest frame [note that 
for self-similar solutions the tetrad-frame connection hr in eq. (pij) is hr = i9t ln(r^'')]: 

QB 

\- Lb — constant . (64) 

6 

Similarly, covariant conservation QnDmT™^ = of dark matter energy implies 
conservation of angular momentum of dark matter in its frame, 

Ld = constant . (65) 

Equation (|65p should be interpreted as meaning that the angular momentum is 
constant along the path of the dark matter, not along the path of the baryons, which 
is different. The dimensionless combination L^/r'^ is however the same in any frame. 

The operator ^'"(9^, equation (j55p . vanishes when acting on any dimensionless 
quantity. Taking times equation (|43aP plus ^* times equation (j436p yields an integral 
of motion for the axial Weyl scalar W, 

^^QB SLb M^ (66) 

Equation (|30cp governing the radial magnetic field scalar B, coupled with 
equation (j46p and the assumption of vanishing conductivity, yields a second order 
ordinary differential equation for B, 

' ' ' i'' dx ' c 



1 

2A 



B + B + 2QW^^ . (67) 



For self-similar solutions, the vierbein coefficient uj is (coordinate and tetrad) 
gauge-invariant, and is related to the axial Weyl scalar W by, from equation ((22l) . 

The vierbein coefficients C™, the tetrad- frame components of the azimuthal Killing 
vector, are not tetrad gauge- invariant. A traditional choice of angular gauge is the 
principal, or Weyl, gauge, defined to be the frame in which the spin-1 components (25) 
of the Weyl tensor vanish. In the principal tetrad, 

where the angular momentum parameter a is 

_ M, dW 

a = — : — . (70) 

GM.r^C^"- dx 

The factor of M,/M, in equation ([69]) brings a to the same scale as the angular 



momentum parameters Ob and Ud of the baryons and dark matter, equation ([38]). The 
factor follows from scaling from conformal time t to Kerr-Newman time iKN at the 
sonic point, where dlnr/dr — 1, whereas 91nr/9<KN = M,/M,, so that 

dT M, , . . 

— at the sonic pomt . (71) 

d^KN M. 

The angular momentum parameter a, of the black hole itself can be defined to be 
the value that would be measured by a distant observer if there were no mass or charge 
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outside the sonic point. Matching the Kerr-Newmann (|496p and self-similar ([66]) values 
of the axial Weyl scalar W at the sonic point yields 

-r'^W + QB QB 

a,M, — = h ifa + id at the sonic point . (72) 

3 6 



4.3. Similarity boundary conditions 

What angular boundary conditions must be adjoined to the boundary conditions of 
the background spherically symmetric solution? 

First, the baryonic and dark matter angular momentum parameters Of, and a^, 
equation l|38[) . must be constants at the sonic point, independent of polar angle 9. 
This means that the angular flow patterns of the infalling baryons and dark matter 
are dipoles, with angular velocity proportional to sin 6*. 

The equation (|67p governing the radial magnetic field scalar i? is a second order 
differential equation, so requires two boundary conditions. However, one boundary 
condition is fixed by the requirement that the solution pass smoothly through the 
horizon, which is a singular point of the differential equation. I use an interative 
shooting method to locate the desired solution. The remaining boundary condition 
on B is in principle set by the requirement that _B, which is proportional to the 
radial component of a dipole magnetic field, should go to zero at infinity (absent 
distant current sources). However, the self-similar solutions generally do not continue 
consistently to infinity. Instead, I simply set B equal to the Kerr-Newman value (|49cp 
at the sonic point. 

The angular momentum parameters Oh and of the baryons and dark matter, 
together with the value of B at the sonic point, determine the total angular momentum 
parameter a, of the black hole through equation (|72p . Since all angular quantities 
scale in proportion to a,, the angular behavior of the solutions is determined in effect 
by a single dimensionless boundary condition, the ratio a^/a;, of the dark matter to 
baryonic angular momentum parameters {B being fixed to its Kerr-Newman value at 
the sonic point). 



4-. 4^. Similarity model 

Figure [1] shows a self-similar solution with accretion rate M, = 0.01, charge-to-mass 
Q»/M, = 0.8, and dark-matter-to-baryon density ratio Pdl Pb = 0.1 at the sonic point. 
This is the same model shown in Figs. 4, 6, and 7 of ^7\. The ratio Od/ab of dark 
matter to baryon angular momentum parameters has been set to 1 far from the black 
hole 

^ = 1 . (73) 
at 

The accretion rate M, = 0.01 (which essentially means that the sonic point 
is expanding at 0.01 of the speed of light, as measured by a distant observer) is 
much larger than any typical astronomical accretion rate. The advantage of a large 
accretion rate is that it is possible to follow the complete evolution of the spacetime 
through inflation to collapse without quantities overflowing numerically. The smaller 
the accretion rate, the more rapidly inflation exponentiates [T7], triggering numerical 
overflow. 

The top panel of Figure [1] shows the various components of the tetrad-frame 
energy-momentum tensor measured in a frame that is center-of-mass in the radial 
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Figure 1. Self-similar model with accretion rate M, = 0.01, charge-to-mass 
Qm/M, = 0.8, and dark-matter-to-baryon density ratio pd/Pb = 0.1, at the 
sonic point. Lines are dashed where a quantity is negative. All quantities are 
in geometric units, c = G = M, = 1. Perturbed quantities (everything except 
T" and T'''") are in units where a, = 1, and should be multiplied by the actual 
a,, a small number. (Top) Components of the tetrad- frame energy-momentum 
tensor T™" in a frame that is center-of-mass in the radial direction and principal 
in the angular direction. (Upper middle) Axial spin-0 Weyl and electromagnetic 
fields W and B, equations 1221 1 and I I28I I. (Lower middle) Vierbein coefficient 
u)M,/M, [the scaling adjusts from conformal time to Kerr-Newman time at the 
sonic point, eq. 1171 l l] . (Bottom) Angular momentum parameters ai, and of the 
baryons and dark matter, equation l|38|l . and a of the black hole in the principal 
frame, equation (I70I I. 
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direction, T**" = 0, and principal in the angular direction, where the spin-1 components 
of the Weyl tensor, equation (25), vanish. The center-of-mass energy density and 
pressure T** and T^^ inflate rapidly at a radius r « 0.4 close to where the inner 
horizon would have been in the absence of inflation. The center-of-mass energy and 
pressure arise almost entirely from the streaming energy and pressure produced by 
the relativistic counter-streaming between the outgoing baryons and the ingoing dark 
matter. The proper energy and pressure of either of the baryon or dark matter streams 
in their own frames change only modestly. 

The angular components T*"^ and T^'^ in the top panel of Figure [T] are plotted 
in units such that the black hole angular parameter is unity, a, = 1, and should be 
multiplied by the true a,, a small number. Thus although the Figure appears to show 
that, in the Kerr-Newman regime r > 0.4, the angular components T**^ and T^^ are 
comparable to the radial components T** and T'"'", in fact the angular components are 
much smaller. In the inflationary regime r < 0.4, the angular components T*''^ and T'''^ 
inflate, but they inflate more slowly than the radial components T" and T^^. Thus 
the angular components, already a small perturbation in the Kerr-Newman regime, 
become an even smaller perturbation in the inflationary regime. 

The upper middle panel of Figure [T] shows the axial Weyl and electromagnetic 
scalars W and B, equations ([22|) and (|28p. while the lower middle panel shows the 
gauge- invariant vierbein perturbation w, equation (|13ap . All three quantities remain 
well-behaved during inflation, showing no sign of any inflationary behavior. This 
demonstrates that the assumption of a slowly-rotating perturbation is consistent: the 
perturbation remains small throughout inflation. 

Note that the uj obtained from equation ([68]) is a conformal angular velocity. 
Figure [T] shows uj multiplied by M,/M,, equation (|7ip . which converts w to a 
true angular velocity that can be compared directly to the Kerr-Newman angular 
velocity (149 a^ . 

The lower panel of Figure [1] shows various angular momentum parameters a. The 
dark matter angular momentum parameter ad is constant, reflecting the fact that the 
dark matter is neutral and pressureless, and therefore falls freely with constant angular 
momentum. The baryonic angular momentum parameter ab decreases inwards, as 
the charged baryons lose angular momentum into the magnetic field, in accordance 
with equation (164]) . The parameter a, equation ([70]) . is defined to be the angular 
momentum parameter in the principal frame. In the Kerr-Newman geometry, this 
angular momentum parameter would be constant, and equal to 1 in the units a, = 1 
used in the Figure. Figure [1] shows that a in the self-similar solution is indeed close 
to 1 in the Kerr-Newman regime r > 0.4. What appears to be a glitch in a at 
r « 0.63 arises because dVF/da; and C, whose ratio yields a, equation ([70|l . go through 
zero at very slightly different radii (in Kerr-Newman, the zeros occur at exactly the 
same radius). No physical divergence is associated with the glitch. All three angular 
momentum parameters remain well-behaved in the inflationary regime: they do not 
inflate. 

The angular structure of the self-similar solution is determined by four equations, 
of which three are algebraic, equations ([64 |) -([66 |) . and the fourth is a second order 
differential equation for the radial magnetic field scalar B, equation (l67|) . Only the 
fourth equation poses a challenge to understanding its behavior. Figure [2] shows that 
in the Kerr-Newman regime r > 0.4, the magnetic B field is closely approximated by 
the Kerr-Newman form (|49cp . Inflation takes place just above the inner horizon, and 
in the inflationary regime r < 0.4, the horizon function A is tiny, and ~ 
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.1 .2 .5 1 
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Figure 2. Radial magnetic field B (left scale), and Regge- Wheeler coordinate 
r* (right scale) for the model shown in Figure [T] The dashed lines superposed on 
the magnetic field B show approximations valid respectively in the Kerr-Newman 
(red) and infiationary (blue) regimes, equations Il49el l and II77I I. 
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B = 



this regime the differential equation ((67|) for B simplifies to 

A_d_ 

^da- 
Introduce the Regge- Wheeler coordinate r* defined by 
d ^ _ A_d_ 
dr * dx 
Then equation (|7^ simplifies to 
d /_d_ 
dr* \dr* 
whose general solution is 

B = Bo + Bie"-"' (77) 

where Bq and Bi are constants. Figure [2] shows a fit to the approximation ([T7| 
with the constants Bq and Bi treated as free parameters, and it is seen that the 
fit is excellent. Analytic and numerical investigation indicates that the constant Bi 
depends on boundary conditions (such as pd/ Pb or ad/ab), and does not tend to some 
universal limit in the limit of small accretion rates (for example, Bi does not vanish). 



2] B^O 



(74) 



(75) 



(76) 



5. Small accretion rate 



The self-similar model of the previous section, 21 has the virtue of yielding exact 
solutions, but it is unrealistic in the sense that it requires a finely-tuned accretion 
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flow, in which the accreting streams have a small dipole angular velocity that remains 
steady over the age of the black hole. 

In fact however, as will now be shown, for sufficiently small accretion rates, 
condition (|8ip , the instantaneous angular motion of the accretion flow has a negligible 
effect on the spacetime geometry of the slowly rotating black hole. For small 
enough accretion rates, the equations describing the angular spacetime structure of 
a slowly-rotating black hole simplify to equations ([82| and ([83| . independent of the 
instantaneous angular motion of the accreting streams. 

For small accretion rates, the interior charge Q can be treated as effectively 
constant as a function of radius, 

Q ~ Q, ~ constant . (78) 

Initially, before inflation, the charge is constant simply because the small accretion 
rate implies a small charge density. During inflation, as first pointed out by [29], 
Lagrangian volume elements remain little distorted even though the tidal force (Weyl 
curvature) exponentiates hugely. The physical reason for the mild distortion is that a 
volume element experiences only a tiny proper time during inflation, so that the tidal 
force, though huge, has too little time to do damage. Because volume elements change 
little during inflation, and because charge density is conserved and the initial charge 
density is small, therefore the charge density remains always small, and the interior 
charge Q is effectively constant. 

For constant interior charge Q, the Einstein equation (j43b[) integrates along the 
worldline of any tetrad frame to 

r^W ^Q.B-Sa.M.+Sn ^ / a{p+p)uW^^ . (79) 

streams "'sonic Pt 

Here the integral with respect to proper time r along the worldline of the tetrad has 
been converted to an integral with respect to radius r using dr/dr = dtr — (3f The 
constant —3a, A/, of integration in the integral ((79)) is established at the sonic point, 
where equation (|72|) holds. 

For small accretion rates, the integral over either of the baryonic or dark matter 
streams in equation ((79|) is small 

a(p + p)MVj«0. (80) 

Pt 

Before inflation, the integral is small because the small accretion rate implies a 
small proper density and pressure p and p. During inflation, the proper density 
and pressure p and p of a stream change only mildly (because Lagrangian volume 
elements distort mildly), and, thanks to conservation of angular momentum, the 
angular momentum parameter a of the stream likewise changes little. The only factor 
that could potentially become large in the integrand of the integral (|80|1 is the ratio 
/ Pt- The integrals in equation ([79]) can be evaluated in any tetrad frame (the same 
frame for both streams) ; one option is to choose the no-going frame, defined to be the 
frame where (3r = 0. Relative to the no-going frame, the streaming velocity vT and 
fit both increase exponentially, but they increase approximately in proportion to each 
other [17]. The result is that the integral ([80]) remains small because p + p \s small, 
and the other factors in the integrand do not spoil that conclusion. 

The integral over either stream in equation (l79|) involves the product a{p + p) of 
the angular momentum parameter and the proper density plus pressure. The integrals 
are small compared to the other terms in equation (|79l) provided that a{p + p) <C a. 
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in geometric units (c ~ G = il/, = 1). If the streams are slowly rotating, a ~ a,, then 
the integrals are small provided that p+p <C 1. But if the accretion rate is sufhciently 
small that it satisfies the stronger condition that 

p + p^a, (geometric units) , (81) 

then the integrals remain small even if the angular momentum parameter a is large, of 
order unity in geometric units. The fact that the integrals remain small, equation (j80p . 
even for large angular momenta means that the instantaneous angular motion of the 
accreting streams has effectively negligble effect on the spacetime structure of the 
slowly-rotating black hole, even in the case that the accreting streams have large 
angular momentum and any angular motion, not just cylindrically symmetric rotation 
about the rotation axis. Of course, ultimately the angular momentum of the black 
hole is equal to the cumulative angular momentum of the accreted streams; but the 
instantaneous angular motion is essentially irrelevant. 

In the approximations ((78| and ((80|) . valid for small enough accretion rates, 
equation ([79|) simplifies to 

r^W w Q.B - 3a.Af. , (82) 

which relates the axial Weyl scalar W to the axial magnetic field scalar B without 
any dependence on the baryon or dark matter streams. 

To complete the equations governing the spacetime structure for small accretion 
rates, equation ()82|1 must be supplemented by an equation for the magnetic field 
B, specifically equation (j30c|) with, since the accretion rate is small, vanishing right 
hand side. Self-similar solutions provide a guide to what happens. As illustrated in 
Figured B is well-approximated by the slowly-rotating Kerr-Newman solution (j49cp 
before inflation sets in, and by approximation (|77p in the inflationary regime, The 
approximation ([77| depends on the Regge- Wheeler coordinate r* , which has the 
property that, as illustrated in Figure [21 the derivative dr*/dlnr goes through 
a maximum near the onset of inflation, and then becomes progressively smaller 
during and after inflation. The net effect is that B flattens out and becomes 
constant during and after inflation, which is true even though the constant Bi in the 
approximation (|77p , which is effectively determined by matching to the Kerr-Newman 
approximation ()49cP at the onset of inflation, is non-zero. For small accretion rates, the 
transition from the Kerr-Newman regime to the regime where B is constant becomes 
increasingly abrupt, and B simplifies to 

(83) 

max(r, r_ j 



where r_ = M, — y/ — Ql is the radius of the erstwhile inner horizon, which is 
destroyed by inflation. The approximation (|83p has been derived on the assumption 
of self-similarity, which is premised on the assumption of constant accretion rate M, . 
However, since the approximation (|83p is independent of the value of the accretion 
rate, it may be guessed that the approximation (|83p remains valid even for a variable 
accretion rate, just so long as the accretion rate remains small. 

The axial spin-0 Weyl scalar is in a fundamental sense the (coordinate and 
tetrad) gauge-invariant quantity that describes the angular part of the gravitational 
field. The Weyl scalar W is related to the angular vierbein coefficients ui and ip by 
equation (|22p . which in the self-similar case becomes equation which depends 

only on lo. Converted to a derivative with respect to Inr instead of x, the self-similar 
equation (|68)) for ui becomes — (/3t/2^'')da;/dlnr = W. In the inflationary regime, the 
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coefficient — /3t/2^'' is so large, for small accretion rates, that dw/dlnr is tiny, so that 
oj is essentially constant as a function of radius r (beware: this statement fails deep 
into the post-inflationary collapse regime). The result is that, for small accretion rates, 
the vierbein coefficient lu is approximately equal to the Kerr-Newman form (j49ap prior 
to inflation, and to a constant during and after inflation: 

^ « f 2A/. . (84) 

max(r, r_)'^ \ max(r, r_)/ 

As with equation ([55]) . equation ([51)1 has been derived on the assumption of self- 
similarity, but in fact the equation is independent of the value of the accretion rate, 
so it may be guessed that the approximation (|84p remains valid even for an accretion 
rate that varies with time, just so long as the accretion rate remains small. 

Equations ((82|) - ([84|) look enticingly simple, but one should be wary of interpreting 
them too literally, especially equation (|84p . For example, if the approximation 
u « constant in the inflationary regime were taken literally, then it would imply, 
from equation ([68]) . that W 0, which is false. Thus it would be incorrect to insert 
the approximation ()84ll into the line-element (I20|) and imagine that the consequent 
Riemann tensor would be valid. The problem is that inflation generates exponentially 
huge numbers, and something that by itself appears sensibly equal to zero may become 
significant when multiplied by a huge number. 

6. Limitation 

The slowly-rotating models considered in this paper have the limitation that the 
rates of accretion of ingoing and outgoing streams are assumed to be constant as a 
function of angular position about the black hole. This limitation is inherited from the 
assumption that the slow rotation is a small perturbation of a background spacetime 
that is spherically symmetric. The previous section showed that the angular motion of 
the accreting streams has negligible effect, for small accretion rates. But the angular 
motion is distinct from the angular dependence of the accretion rate. 

It would be possible to consider, in addition to the slowly rotating perturbation 
considered in this paper, perturbations to the accretion rate as a function of angular 
position over the black hole. This could be the subject of future investigation. 

7. Summary 

The internal structure of a slowly rotating black hole that is undergoing inflation at 
(just above) its inner horizon has been derived. The approach has been to introduce 
a slowly-rotating perturbation to a spherically symmetric black hole undergoing 
inflation. It has been shown that the rotational perturbation remains small throughout 
inflation, so the perturbation assumption is self-consistent. 

The equations governing the angular behavior, Sj^l decouple from the radial 
behavior. Consequently all of the conclusions concerning inflation in spherically 
symmetric black holes carry through unchanged for slowly rotating black holes. This 
conclusion supports [Tj's conjecture that inflation in rotating black holes is probably 
similar to inflation in spherical black holes. 

Exact self-similar solutions for slowly rotating black holes have been obtained, 
21 The self-similar solutions require a special accretion flow, in which the accreting 
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streams have a small dipole angular velocity that remains steady over the age of the 
black hole. 

It has been shown, SjSl that for sufficiently small accretion rates the instantaneous 
angular motion of the accretion flow has a negligible effect on the spacetime geometry 
of the slowly rotating black hole. The angular structure of the black hole is described 
by the axial spin-0 Weyl and magnetic field scalars W and B, which for small 
accretion rates are given approximately by equations (|82p and (|83p . independent of 
the instantaneous angular motion of the accretion flow. 

The slowly-rotating solutions obtained in this paper all have an accretion rate 
that is constant as a function of angular position over the black hole, a limitation 
inherited from the assumption that the slow rotation is a small perturbation from 
spherical symmetry. 
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Appendix A. Finite conductivity 



This Appendix gives, for reference, the self-similar equations for the case of non- 
vanishing baryonic electrical conductivity a. Equation (jA.3[) replaces equation (j64[) . 
while equation (jA.4[) replaces equation (|67p . 

To admit self-similar behavior, the conductivity must have dimension a cx r^^. 
If the conductivity is assumed to depend only on the baryon proper density pb, then, 
since pb oc r~^, it follows that [the following repeats eq. (48) of [T5] ] 



Kp 



1/2 



" - (4^)1/2 ' (^-^^ 

where k is a phenomenological dimensionless conductivity coefficient. A dimensionless 
conductivity s can be defined as [the following repeats eq. (51) of [TS]] 

s = Anar = K.{A'Kr'^ pf/'^ . (A.2) 

In terms of the dimensionless electrical conductivity s, the equation expressing 
conservation of the angular momentum of baryons is 

A (OR , r\ . -f^ (A.3) 



Ax \ & J 3 
and the equation governing the radial magnetic field scalar B is 



1 

2A 



-1 



A d 

^dx 



B + B + 2QW 



IQ_ 

rSj' 



(l + ^e*) 



sdB 
2~dx 



(A.4) 
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